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We present explicit expressions for a large set of hierarchy wave functions on the torus. Included 
are the Laughlin states, the states in the positive Jain series, and recently observed states at e.g. 
v = 4/11. The techniques we use constitute a nontrivial extension of the conformal field theory 
methods developed earlier to construct the corresponding wave functions in disc geometry. 
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I. INTRODUCTION 

A basic concept in the physics of the fractional quan- 
tum Hall effect (FQHE) is that of an incompressible liq- 
uid. Theoretically, there are two main types of such liq- 
uids depending on the quantum statistics of the quasi- 
particles: the abelian and the non-abelian liquids. Of 
these, only the former are firmly established experimen- 
tally. The abelian liquids are further divided into single- 
and multi-component ones, where the former include the 
prominent filling fractions corresponding to the Laughlin 
states at v = l/(2q + l)[l[ and the composite fermion 
states in the Jain series v = n/(2qn± 1)@. Multicompo- 
nent liquids can form in systems with spin or pseudospin 
degrees of freedom, i.e. in partially polarized or bilayer 
systems Q. 

There has been two major approaches to the abelian 
single-component quantum Hall (QH) liquids, the 
Haldane-Halperin hierarchy 0, @L an d the composite 
fermion (CF) approach of Jain [2, [||. In the former 
scheme, the quasi-particles that exist in the vicinity of 
a given Laughlin state condense into a new Laughlin- 
like state; this procedure can be repeated thus forming 
an ever more complicated hierarchy of coexisting incom- 
pressible liquids, while in the CF scheme the prominent 
fractions are understood in terms of filled effective CF 
Landau levels. A major advantage of the CF approach 
is that it provides explicit, and numerically very accu- 
rate, wave functions, while an appealing feature of the 
hierarchy scheme is that it treats all filling fractions on 
the same footing, including the newly observed ones Q 
at e.g. v = 4/11. 

In a series of recent papers [1, H, Ell EL E2| we have 
given a concrete realization of the Haldane-Halperin hi- 
erarchy both by providing a large set of explicit and 
testable wave functions derived using conformal field the- 
ory (CFT) techniques, and by the construction of an ex- 
actly solvable model describing interacting electrons in 
the lowest Landau level on a thin torus. The CFT con- 
struction yields candidate wave functions for every filling 
fraction obtained by successive condensation of quasielec- 



trons (as opposed to quasiholes) that reduce to the ex- 
act solutions on the thin torus and coincide with those 
due to Laughlin and Jain, whenever these exists. The 
Laughlin states are the exact ground states of a short 
range model Hamiltonian, but no such Hamiltonian is 
known for any other hierarchical state. Both the Laugh- 
lin and Jain wave functions are however known to be ex- 
cellent approximations to the numerically found Coulomb 
groundstates for small systems. 

The standard way to quantitatively decide if a pro- 
posed QH wave function captures the physics of a given 
phase is to study small systems numerically, typically 
by exact diagonalization, and compare the numerically 
obtained energies and wave functions with the proposed 
ones. To carry out such comparisons one has to choose a 
geometry. Early calculations were performed in a physi- 
cally motivated disc geometry, but rather than using an 
explicit confining potential, the total angular momentum 
was fixed. The drawback of this method is that small sys- 
tems might exhibit large, and unphysical, edge effects. 
These can be eliminated by using a finite geometry, and 
most large scale calculations have been performed on the 
sphere. Wave functions on the torus, on the other hand, 
turn out to be more complicated, but this geometry pro- 
vides other advantages. Not only are there no edges, 
but also no curvature, and since the torus is a genus 
one manifold, the topological properties of the different 
phases will be reflected in the ground state degeneracy 
[l3j . Also, on a torus the one-dimensional nature of a 
Landau level is explicit: there is a natural mapping of the 
two-dimensional problem onto a one-dimensional lattice 
model that is solvable on the thin torus. There has been 
considerable progress in understanding various phases of 
the QH system in terms of exactly solvable models that 
describe the thin torus [1, [Q, Oj, OJ, OJ, El, for 
details see Ref. |l9|. 

Altogether we think that there are strong reasons to 
construct the hierarchy wave functions on the torus: it 
will facilitate comparisons with numerical calculations, 
via the ground state degeneracies it will give information 
about the topological order, and it will provide a way to 
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extend the exact results in the solvable thin torus limit 
to the experimentally accessible regime corresponding to 
a thick torus. Explicit wave functions on the torus have 
previously only been obtained in a few special cases where 
the structure is simple (20l . [2lj , and in particular, the 
torus version of the Jain, or the hierarchy wave functions 
have, to our knowledge, not been constructed. 

In the present article we extend the CFT construction 
to obtain explicit wave functions in the torus geometry. 
In Section [IT] we briefly review the necessary CFT ma- 
chinery, as well as some technical details about the low- 
est Landau level on a torus. In Section 1X111 we use this to 
construct wave functions on the torus. We analyze the 
Laughlin states (one field) and the v = 2/5 state (two 
fields) in some detail, and then provide the basic relations 
and logical steps needed to derive the wave function for 
a general hierarchy state that is obtained form a Laugh- 
lin state by successive condensations of quasielectrons. In 
Section[V]we describe an alternative choice of charge vec- 
tors that simplifies the structure of the wave functions, 
and is very useful for explicit calculations, although it 
is not suitable for deriving the general formulas of the 
previous section. In Section [VII we compare the v = 2/5 
wave function with the result of an exact numerical diag- 
onalization by calculating overlaps. Section IVIII includes 
some concluding remarks and an outlook. Several techni- 
cal steps are explained in Appendices O and [B] and as a 
service to the practically oriented reader we also give ex- 
plicit expressions for the states at v = 3/7 and v = 4/11 
in Appendix [C] 

II. HIERARCHY STATES AS CFT 
CORRELATORS 

In this section we review the procedure, for express- 
ing QH wave functions in a disk geometry as anti- 
symmetrized sums of CFT correlators, and emphasize 
several points that will be important for the generaliza- 
tion to the torus. 



A. General considerations and the Laughlin state 

Our aim is to relate the hierarchy wave functions on 
the torus to the holomorphic conformal blocks that build 
the correlation functions of vertex operators in certain 
CFT's. We first briefly recapitulate the construction in 
the plane outlined in Refs. 0, [ll|. 

The relevant CFT's have actions 

S[<p a ] = £ Jd 2 xd^ a d^ a , (1) 

where fi = 1,2, the metric is Euclidean, g is an over- 
all normalization to be fixed later, and the boson fields, 
ip a (z,z), with z = x\ + ix2 — x + iy, ar e compactified 
on circles with radii R a . In loose analogy to the stan- 
dard CF construction we will need n fields to describe 



the CF state with n (partially) filled CF-Landau levels 
(27| . Note that the compactification implies that tp a are 
angular variables, so on manifolds of higher genus, there 
are field configurations with nontrivial winding around 
the different handles. 

The primary fields are given by the vertex operators, 

V Q {z,z) =: e ^=i<2"<M^) : , (2) 

where the colons, which we suppress in the following, de- 
note normal ordering, and Q is an n-dimensional charge 
vector, 

where are integers. 

On the plane we choose the normalization g in the 
action so that the two-point function of the scalar fields 
is given by 

(ip a (z,z)ip b (w,w)) = -S ab \n\z - w\ 2 , (4) 

and hence a contraction of two vertex operators is given 
by 

(V QM (z,z)V Q(j) (w,w)) = |z- W | 2Q<!) - Q<J \ (5) 

for any charge vectors QW and Q^. 

The statistics of the particles is coded in the opera- 
tor product expansion (OPE), or equivalently, in the two 
point function, of the corresponding holomorphic vertex 
operator Vq(z). Setting i = i in ([5]). we find that the 
holomorphic part of the vertex operators with |QW | 2 odd 
has fermionic statistics. On the torus, the analysis is 
more involved, but the local properties of the vertex op- 
erators are the same as on the plane as is evident from the 
pertinent OPE's. The statistics of the quasiholes [22j and 
quasielectrons Q can be understood in a similar manner. 

The QH wave functions can be obtained as correlators 
of the holomorphic part of the vertex operators. The 
simplest example is that of a Laughlin state at v = l/q. 
Here we have only one field, ip, with radius R 2 = q and 
charge e = q, and get 

N 

2—1 i<j 

= * Tpy q (zi), (6) 

where Ob g is a constant neutralizing background, and $ 
a singular phase factor that can be properly defined by 
replacing the continuous background charge by a lattice 
of thin flux tubes. Then $ = $(0j(n o )), where <f>i(n a ) 
is the relative angle between the coordinate zi and the 
lattice vector n a . For the details of this procedure, which 
can be generalized to the full Jain series, see Appendix 
A in Ref. ©. 
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Alternatively, we can recover the wave functions by 
factoring the full correlation function into a holomor- 
phic and antiholomorphic part, also called conformal 
blocks, each accompanied by the square root of the non- 
factorizable exponential factor 

TV 

tfl V^izu - Zi )O hg ) = n |zy \ 2 "e-^ ^ ^ 



1=1 



1<J 



and are constructed recursively by 



Vi 



V a+ i = d z V a e 



1, 



,,n- 



1(8) 



where 71 = V2&i + 1, 7«+i = V 2fc Q +i - 7q and <A*+i 
is a new bosonic field. The electrons are divided into n 
sets I a of size M a and the wave function is given by 



n 



z « e -;d*£il*l a 



n 



(^l/ 9 (2!i))^l/ ? (Zi) 



(7) 



where = £j — z^. It is this latter procedure that will 
generalize to the torus. The background charge must be 
included to obtain a non-vanishing correlator, and will 
be treated in detail in Section ITTT1 For now it suffices to 
mention that it gives rise to the correct nonholomorphic 
dependence, and nothing more, as it does in the disc 
geometry. 

On the torus, correlation functions will not factorize as 
111 0, but are given by an (in general infinite) sum over 
such terms. For a special class of CFT's, this infinite sum 
can be rewritten as a finite sum over extended conformal 
blocks. Such theories are called rational CFT's [23j. For 
a bosonic action of the type (JTJ) this requires that the 
radii, R a , are of the form = 2p/p' with p and p' 
relatively prime. This particular rational CFT is called 
a rational torus. Clearly the Laughlin state considered 
above falls into this class, as do the hierarchy states that 
are discussed in the following section. 

In the general case we will thus extract not only a sin- 
gle wave function, but a whole set. This provides an 
important consistency check on our method, since the 
degeneracy of a lowest Landau level (LLL) state on the 
torus is known from the general symmetry analysis due 
to Haldane (24[. The degeneracy of a state with filling 
fraction v = p/q, includes a factor four related to choos- 
ing periodic or antiperiodic boundary conditions along 
the cycles of the torus, and a g-fold degeneracy related 
to the position of the center-of-mass. As we will see, our 
construction precisely recovers this degeneracy. 



B. The hierarchy wave functions 

The hierarchy wave functions at level n, obtained by 
successive condensation of quasiparticles only, are de- 
scribed by a set of integers {ki, . . . , fc n }, related to the 
densities of the quasiparticle condensates. In the CFT 
scheme of Ref. 11 these wave functions are constructed 
from correlators involving n different vertex operators, 
V a ■ They depend on the boson fields (f a , a = 1 , . . . , a 



(9) 



where A denotes antisymmetrisation and (...) the cor- 
relation function in a suitable background field. The M a 
are determined recursively as explained in Ref. [l2l . The 
filling factor of '5, v n — p n /q n , is obtained recursively 
from 



Pn 2k n p 

n—1 Pn—2 

q n 2k n q n -i — q n -2 



(10) 



with the initial conditions qq = pa = 1, and the Laugh- 
lin series given by v\ = l/(2fci + 1), ki integer. Taking 
kj = 1 amounts to having maximal density in the j th 
condensate, and the Jain series v n = n/{2nki + 1) corre- 
sponds to choosing ki = • • • = k n = 1. 

To evaluate the correlators of the operators ([8]) on the 
torus, we will use a different representation of the vertex 
operators, where the charge vectors are explicit: 



V, 



(11) 



Here, are derivative operators to be discussed be- 

low. A vertex operator without a hat is a descendant, 
whereas a hat marks a primary field. The charge vec- 
tors Q( Q '") are not uniquely defined. We will mostly 
use a parametrization where the hierarchy construction 
is manifest: 

Q («,n) = ( j* ...,|2=i,|L ,0, ...,0). (12) 
til tl a -i K a 

The CFT charges, Cj and qj, and the compactification 
radii, Rj, are connected to the denominators of the filling 
fractions by 



Qj - Qj-i 



Ri 



(13) 



For the Jain series, for example, all charges are equal to 
2k\. When there is no ambiguity, we suppress the level 
index n on the charge vector and write just Q^"'. 

Other choices of basis sets represent the FQH state 
equally well as long as the bosonic fields, ipi,...,ip n , have 
rational R? and the inner product of two vectors is un- 
changed: Q( Q ) • Q( Q ') = <W + 2h + YT,=2 2 (% " 1)> for 
a < a' . In the special case of the Jain series, this rela- 
tion simplifies to: ■ Q( q ') = 6 aa > + 2k x [|. This is 
discussed in more detail in Section fVl 
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For later computations in Section IIV C[ we will also 
need the following relation between the number of states, 
N a , and the charge vectors: 



Va : ^Q (Q) ■ QWMp = N s 
0=i 



(14) 



Together with "Y^pMp = N = vN s , this gives a consis- 
tency condition on N s . 



III. QH STATES ON THE TORUS - GENERAL 
CONSIDERATIONS 



In this section we describe the techniques we use to 
construct the hierarchy wave functions on the torus. 



A. Lowest Landau level on the torus 

We consider a system of N charges — e on the torus 
with periods L\ and L 2 . In Landau gauge, a homoge- 
neous external magnetic field perpendicular to the sur- 
face is described in terms of a vector potential A = Byk. 
Most of the ti me we w ill set the magnetic length equal 
to unity, I — ^Jhc/eB = 1 

In the presence of an external magnetic field, consistent 
periodic boundary conditions can only be enforced up to 
a gauge transformation. A natural way to do this is to 
use the magnetic translation operators [HI 



t(I) = eHV-iAj-iixx ) 



(15) 



where 1 parametrizes the translation, and x the electron 
coordinate. In Landau gauge we define 



t 



t a = t(fty) = e 



(16) 



where N s = LiL 2 /2ir e Z is the number of flux quanta 
penetrating the surface of the torus. The translations 



and t 2 s commute. Hence, periodic boundary conditions 
can be consistently formulated as 



exp (i<j> M ) ip ; [1 = 1,2 



(17) 



where <j)^ can be interpreted as solenoid fluxes through 
the handles of the torus. Once the fluxes are fixed, 
the same boundary conditions apply for all states in the 
Hilbert space of the system. 

A complete set of commuting operators for a v = 
p/q QH system, with p and q relative primes, is given 

by {HfTifT^}, where T M = n»=i*/M are centcr-of- 
mass translations and if is a translationally invariant 
electron-electron interaction. The states are labeled by 
{E, Ki,K%} where the quantum numbers K a of the sys- 
tem relate to the eigenvalues of T\ and T%, which are 
given by e 2viKa ^ Ns respectively. 



The lowest Landau level wave functions, in Landau 
gauge, on the torus have m any special properties dis- 
cussed for instance in Refs. [2J, [25|. They all factorize 
as 



,z N ) 



f(zi, 



z N )e 



(18) 



where / is holomorphic in all its arguments and is given 
in terms of generalized -^-functions 



(*|r)= J2 e 



1-77 T (fe+a) 2 2-rri(k+a)(z+b) 



(19) 



k=-c 



For example, a Jastrow-type factor rii<j — z j) q on ^ ne 
plane becomes Y\ i< j°^i{- i j^ 1 -\ij^) q on the torus, where 
$1 corresponds to a = b = 1/2. 

From the boundary conditions in Eq. \17\ one derives 
the corresponding conditions on the holomorphic func- 
tions / as will be discussed in detail below. 

Because of translational invariance, the center-of-mass 
(CM) dependence of the wave function can be separated: 



f(zi, ...,z N ) = F cm (Z)f r& i{z\, . . .,zn), 



(20) 



where Z = J^. Zi/ L\ and f re i is independent of Z . While 
for the Laughlin states f re i is uniquely determined, (at 
least on the plane, sphere and torus), by its leading short- 
distance behaviour, this is not the case for the Jain states, 
or more general hierarchy states. Of course, once the 
relative part is fixed, the degeneracy of the corresponding 
level is given by the number of center-of-mass functions 
compatible with the boundary conditions. 

For hierarchy wave functions, we have not been able to 
explicitly separate the CM part even in the simple case 
of the Jain sequence. Here, the CFT techniques have 
proved to be extremely useful in that they allow for a 
direct construction of the wave functions without ever 
separating the CM part. 



B. The background charge 

A careful calculation shows that the correlators 
Vq(t) {zi, Zi)) vanish unless they satisfy the charge- 
neutrality condition J2i Q^" 1 = 0- O n the plane, there 
are two standard choices for the neutralizing background. 
For N particles with charge vectors, Q a , one can ei- 
ther assume a compensating charge at infinity, Ob g = 



~ iN J2 a Qa<Pa(.Zao,Zao) 



or a homogeneous droplet, 



O 



by 



JJ exp[-zi? a p a / d 2 x(p a (z, z)\ , (21) 



with p a = po/Ri, Po = 1/(2tt£ 2 ), being the density of 
the o th set. The latter method, which is more physical, 
is the one used in ([7]); it reproduces the correct non- 
holomorphic dependence of the LLL wave functions in 
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the symmetric gauge. On the torus, the first choice is 
not possible at all, therefore we will use the latter. 

As in the plane, this homogeneous distribution cannot 
be realized using the operator content of the bosonic ra- 
tional CFT, which implies that, in general, we cannot 
expect the conformal correlation functions to factorize 
into chiral components. On the plane, this did not cause 
any problems. On the contrary, the homogeneous back- 
ground charge actually contributed the gaussian non- 
holomorphic dependence needed in the LLL wave func- 
tions. Below, we find the same result on the torus. 



C. The derivative operators X> ( ' 

One of the basic difficulties in translating the hierar- 
chy wave functions from the plane to the torus is related 
to the meaning of the derivative operators T>( k \ In the 
plane, they are the holomorphic partial derivatives d\ 
implying that the operators V a are descendants of the 
primary fields in the theory. Using standard techniques, 
the corresponding correlation functions can be written as 
a product of derivatives acting on the correlation func- 
tion of primary fields only. As discussed in Ref. the 
derivatives act only on the polynomial part of the wave 
function, that is on the conformal block, and not on the 
exponential factor. Thus, in the plane, the calculation of 
the wave functions reduces to the calculation of correla- 
tors of primary fields. 

The origin of the derivative operators is easily under- 
stood in terms of composite fcrmions, where they are the 
remnants of the effective Landau level structure. Pro- 
jection onto the lowest Landau level converts the z^'s, 
which occur in higher Landau levels, into derivatives, d%. 
In particular, this means that if a derivative is excluded, 
the wave function will vanish since it amounts to putting 
a particle in an already filled Landau level. For the gen- 
eral hierarchy state, the situation is more complicated, 
but numerical calculations indicate that the derivatives 
are necessary also here. 

Turning to the torus, we have not been able to prove 
that the derivatives are needed in order to get non- 
vanishing wave functions. On the other hand, this is al- 
most obvious since for large tori, the wave functions must 
be very similar to those on the plane, and these do need 
derivatives. We have also checked that the correlators of 
the primary fields do vanish after antisymmctrization in 
simple cases, e.g. for v = 2/5. 

A natural way to construct our wave functions on the 
torus, would thus be to supplement an appropriately pe- 
riodized version of the Jastrow factors occurring in the 
disc version, with derivatives and a suitably chosen CM 
part. Such an ansatz, however, does not work. A deriva- 
tive acting only on the holomorphic part of a correlator 
destroys the quasi-periodicity of the wave function. [28| 

We conclude that we must construct the T>( k > 's in (ITTj) 
from operators that preserve the Landau level index and 
the boundary conditions, ie they must commute with the 



one-particle hamiltonian and the magnetic translations, 
tjli, around the cycles of the torus. Requiring, in ad- 
dition, that they approach ad + bd + c (where a, b, c are 
constants) on a large torus, leads to the finite translations 
in the x-direction: 

t% = = e^ 9 ^ , (22) 

where there are at most N s distinct possibilities for the 
integer k. We note that these operators will preserve the 
quantum numbers of the many-body wave functions. A 
more thorough analysis shows that, in fact, only \N s /2 — 
1] — that is, N s /2 — 1 if iV2 is an even number and (N s — 
l)/2 if N s is odd — are independent. Forming from 
linear combinations of the operators (|2"2"j) . 

\N s /2-t] 

v {k) = J2 ( 23 ) 

771=0 

we can first compute the correlators of the primary fields, 
and then act with the 2?( fe )'s, just as in the plane. To get 
the correct limit for a large torus, it is natural to demand 
that X>W d k as N s = L 1 L 2 /{2Tri 2 ) -> oo, but this is 

(k) 

not enough to determine the coefficients c m uniquely, 
and may in fact not even be necessary. [29( For the states 
at level two and three that we have tested, we used, 

X> (1) = t x - 1 (24) 
X> (2) = t\ - 2t x + 1 . 

(The constant in T>^ can be ignored as it gives no con- 

(k) 

tribution to wave functions.) That the Cm s cannot be 
determined uniquely might seem a serious drawback, but 
in Section IVII where we present numerical tests of our 
wave funtions, we will see that for the v = 2/5 state, the 
simplest choice, X>W = t±, already gives a good result. 
Furthermore, we show how this can be improved to yield 
a very good wave function by taking a linear combina- 
tion of states where Z>W = t\, k small. We believe that 
adding such contributions from higher values of k is the 
torus counterpart of correcting the wave functions in the 
plane by adding contributions involving descendants of 
the primary fields defining the representative wave func- 
tions 1221. 



D. Expressions for the correlators 

Following the strategy outlined above we now ignore 
the derivative operator in V a , and calculate the correla- 
tors of the corresponding primary operators, V a , in the 
presence of a homogeneous background. Since the scalars 
(p a are not coupled, the correlation function © factor- 
izes and can be calculated using a straightforward gen- 
eralization of known techniques (see e.g. Chapter 12 of 
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Ref. \2 



n n vqwte^jOfc) 

a=l i a el a 
n 

< n n e^^^Mo bg ) 

a=l i a el a 
n n 

n< n n e ^w^j 0b9) 

a— 1 a— 1 i a Gl a 
n 

\^as{z l] )\ 2 \{^\Z^\Z^)f^\z ll z l ) (25) 



with 



vw^o = n n ^(^w 9W - qw 

7=1 i^<j~ f el~ f 

n 

x n n *i(^w q(b) - qw . (jw) 

a</3 



^(0)^(0)^(0)) = e -2«p a fl a /^ [a e , mZ -a e , m 2] 



e.m— — 00 



x e M-< m --< J e H Q - m ^- a =''" ^r] (27) 

and 

/ff(^) (28) 

= n n e~ paRaQ(cr) 1 d2x ^w^^-^-w , 

a = l i a £/ a 

where 21^3 = (zj a — zjJ)/Li etc., and 



a=l 



(29) 



a— o+l 



with Z Q = ^ gJ z ia /Li being the (dimensionless) CM 
coordinate of the electrons in the set 7 Q . We also intro- 
duced the notation r = iL2/L\ for the modular param- 
eter describing the torus, and a e ^ m = e/R a + mR a /2, 
dt e ,m = e/R a — mR a /2, where e, m are integers, to 
parametrize the electric and magnetic sectors of the CFT 
Hilbert space. Note that depends on the filling frac- 
tion v n . In order to avoid confusion, we will often indicate 
the level of hierarchy by the subscript n, Zf?' = Z^ a \ 
When an explicit expression is needed we will instead 
denote it by the filling fraction, like e.g. Z^X. 

The first exponential in (|27p comes from the effect of 
the background charge for the different fields and van- 
ishes, if the integration domain is chosen to be J d 2 x = 



f-Li/2 ^ x S-L2/2 ^ e m teg ra l m can be calcu- 
lated exactly, and using the same integration domain, we 
obtain, 



In 



Li 



7(0,0) + 27r?/ 2 .(30) 



In appendix A, we show that a different choice of in- 
tegration domain only amounts to a coordinate shift in 
the wave function. This, however, is not obvious, but a 
rather non-trivial consequence of the homogeneity of the 
states. 

Substituting (|30|) in ([28| yields, up to a constant fac- 
tor, 



/£>(*,*) =n n 



-~2irR aPa 



(31) 



a=l i a el a 



Finally we can use the relation Y^a=i I Ra = 1, which 
simply expresses that all electrons have unit charge [9|], 
to get 



(32) 



which is exactly the non-holomorphic gaussian factor ap- 
propriate to Landau gauge. 

Since the background charge does not alter the form of 
the charge-lattice summation the conformal blocks have 
the same structure as in a rational CFT. We stress that 
because of the gaussian factor fb g , and that alone, the 
terms in the sum that give the full correlator, cannot be 
factorized into holomorphic and antiholomorphic blocks 
as it would should the operator content be purely that of 
a rational CFT. As in J7J), we extract the QH wave func- 
tion as the holomorphic conformal block times the square 
root of the gaussian factor. In the following, we assume 
a homogeneous neutralizing background and only write 
the fully (anti)holomorphic blocks omitting the gaussian 
factor. 



E. Charge lattice sums 

In order to diagonalize magnetic translations in the 
Hilbert space spanned by the holomorphic conformal 
blocks, it is useful to simplify the charge-lattice sum. If 
the compactification radius is of the form R 2 — 2p/p', 
then we can express the infinite sum in (|27|) as a finite 
sum over extended conformal blocks (see Appendix IB]) 



E 



p'-l 2p-l 

E FrA zia) )F-rAZ [a) ) 1 (33) 

r=0 s=0 
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where the 2pp' functions 



(2pp' k-\-rp-\-sp') 2 /2pp' 



X e p 



(2pp' ' k-\~rp-\- sp' ') Z^ 



(34) 



span the Hilbert space of center-of-mass motion. Under 
single-particle lattice translations they transform as fol- 
lows: 

F r,s{Z^ + c) _ 2-kics/R 2 
Fr,s + c(Z^) 



(35) 



Here c are the integers related to the charges in the charge 
vectors 



IV. 



QH STATES ON THE TORUS - EXPLICIT 
CONSTRUCTIONS 



We now have all the pieces needed to extract the wave 
functions — except for one essential ingredient: if we ex- 
tract the basis functions ip s from the correlators using the 
recipe J2 s [4's(z i )]*'t(; s (zi) discussed in Section [U these 
will not satisfy the boundary conditions (fT7|) . In this 
section, we describe in detail how to construct the cor- 
rect linear combinations for the hierarchy states at filling 
fractions v n (fTT)|) . 

We start by constructing the wave functions explic- 
itly in two simple cases: First the Laughlin state at 
v = l/(2k + 1) and then the simplest level two state, 
the v = 2/5 in the Jain sequence. The reader who care- 
fully studies these examples should get a pretty good idea 
about the strategy for attacking the general case, the de- 
tails of which are given in the last subsection. 



A. The Laughlin state 

As already discussed, the Laughlin wave function for 
v = l/q on the torus [20| can be extracted from the Ap- 
point correlation function of the vertex operators Vq(i) 
depending on a single field (pi compactified on a circle 
with radius R\ — q (i.e. p — q and p' — 2), and with 
charge vector = q/R\ = ^fq. In this case there is 
only one condensate and = 3Z. 

Using (|23|) and (|3"3")) we can express the holomorphic 
wave functions in a basis given by the chiral conformal 
blocks, 



A >S (*) = I] *1 l T )" e " 1 ^r,s (Z W ) 



(36) 



i<3 



A maximal linearly independent set of states is obtained 
with r = 0, 1 and s = 0, 1, . . . , 2q— 1. The number of con- 
formal blocks is 2pp' — Aq which is the expected number 



of degenerate states — a factor q coming from the CM 
degeneracy, and a factor 2x2 from the different bound- 
ary conditions (fl7| . The q-fold degenerate multiplet of 
Laughlin wave functions is given by the linear combina- 
tions of the basis states, (|36|) . that transform irreducibly 
under magnetic lattice translations and satisfy the same 
solenoid flux conditions. 

To find the physical states we have to diagonalize the 
action of the single particle magnetic translation opera- 
tors, t£ j, on the full wave functions. To make the anal- 
ysis more transparent, however, let us first summarize 
the relevant transformation properties of the holomor- 
phic components of the basis states. The odd Jacobi 
theta function transforms as 



Mz + 1\t) = -0i(z|t) 
0i(« + t|t) = -e-' l7rT e- 2mz 



(37) 



The transformation properties of the holomorphic center- 
of-mass functions can be read off from (|35|) and are given 
by 

F r , s (Z (1) +q) = (-l) r FrAZ (1) ) 



TrAZ {1) + qr) = e 



— inrq — 2-niZ^ 



F r , s+q (zW); (38) 



the functions should also satisfy the periodicity condi- 
tions 



T r , s+2q {Z^) = F r ,s{Z {1) ), 
T r+i AZ {1) ) - TrAZ {1) )- 



(39) 



From these relations, we can work out the action of the 
magnetic translations on the basis functions (|3"6"|) ; in the 
x-direction this gives: 



n 



+ l\T) q TrAZ (1) + q) 



-1) 



q(N-l)+r 



(40) 



and in the y-direction: 

^2,1 Vv,s ( z i) 

yr djjzij + T\T) q T r AZ (1) + qr) ;lmr 2l + ±(ImT) 2 



(-1) 



q(N-l) 



, (41) 



where, in the last step, we used the identities Nq = N s 



and LiL 2 — 2nN s . This shows that ip r s are eigen- 
states with eigenvalues independent of s, while t 2 f maps 
the function ip r}S (zi) into the linearly independent func- 
tion ?p r ^ s+q (zi). Since, however, both transformations are 
independent of s, we can, in a unique way, satisfy the 
boundary conditions (| L7[) by forming the linear combina- 
tions 

H%(ZV) = T r AZ {l) ) + {-iyT rrs+q (Z^) (42) 



= EC- 1 : 



kt iTTTq(k-\-ai ) 2 2iri(k-\-ai)qZ 
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where ai = - + £, t = 0, 1 and s = 0, . . . , q — 1, which 
amounts to a change of basis for the conformal blocks 
spanning the Hilbert space of the CM wave function. 
Thus, defining 



1>r t t,te) =]jMzij\T) q Hi%(ZW)e-2^ (43) 



we finally have 

^■VV,m(^) = (-l) q{N - 1)+r ^r,t,s{Zi) 
$&r,t,te) = (-l) 9(JV - 1)+ Vr,t, S -(^) ■ (44) 

Thus, the states ip r ,t,s are eigenstates of t^\, a — 1,2. 
For a fixed number of particles, the four different choices 
of the solenoid fluxes in Eq. [T7] precisely correspond to 
the four combinations of the quantum numbers r and t, 
while the quantum number s determines the position of 
the CM. For a translationally invariant hamiltonian, this 
implies a g-fold degeneracy. 

The many-body quantum numbers are related to the 



CM translations T x = U.f=ih,i and T q = Il*Li*2,i 
which, together with the hamiltonian, form a maximal 
set of commuting operators. The operator T% acts only 
on the holomorphic CM piece, and we get 



riV 



= {-l) r e 2 ™' q H { rl E {Z {1) 



(45) 



Comparing with the definition of the quantum numbers 
K h given by = e 2wiK '/ N ^ip, it follows that K x = 
(rN s /2 + Ns) mod N s . T 2 acts on both the gaussian 
and on the holomorphic CM piece. We get 



n 



H?1- S {ZV 
= (-1) 



H^JZW), (46) 



and hence K2 = (tN s /2) mod N s . This analysis estab- 
lishes that states with unequal s have different quantum 
numbers and are therefore orthogonal. We conclude, in 
agreement with Ref. HO, that the degeneracy of Laugh- 
lin's states, as obtained from conformal correlators, is 
given by the denominator of the filling factor. 

Let us summarize the result for the case cf> a = 0, where 
r = t = (N s — q) mod 2. In order to write the result in 
a more conventional notation, we note that the holomor- 
phic CM function can be written in terms of the Jacobi 
theta function (fl"9)) . Comparison with (|42|) gives 







(N a -q)/2+ s/q 
(N s - q)/2 



(qZ\qr). (47) 
In this notation, the g-fold degenerate multiplet of 



Laughlin wave functions is given by 

JV 



til 



X d 



(N s -q)/2 + s/q 
(N s - q)/2 



(qZ\q T ), (48) 



with s = 0, . . . , q— 1. Up to an overall constant the result 
is identical to that obtained in Ref. [2(1 



B. The v = 2/5 state 

As a first example of a second level hierarchy state, we 
now explicitly compute the wave function for the filling 
fraction v = 2/5. We use the same approach as for the 
Laughlin wave functions: the holomorphic part of the 
correlator, together with the square root of the gaussian, 
yields a set of basis states. We construct all linear com- 
binations which are eigenfunctions of the single-particle 
magnetic translation operators and t 2 1 . For simplic- 
ity, we set both solenoid fluxes to zero. This yields five 
candidate wave functions with T\ eigenvalues e 2mnv . 

The correlator is built from vertex operators con- 
structed from two bosonic fields with radii i?f = 3 and 
i?| = 15. Since this is a Jain state, corresponding to 
two completely filled CF Landau levels, the sets contain 
equal number of electrons so that M\ = M2 = N/2. In 
the hierarchy picture this amounts to having a maximally 
dense condensate of quasielectrons on top of the v = 1/3 
parent state. 

The charge vectors are given by 



Q (2) = (4= ' 



(49) 



V^' a/15 



and we get Z^ = 3Zi + 2Z 2 and Z^ = 5Z 2 where Z x 
and Z2 are the CM coordinates of the two sets respec- 
tively. Following the strategy outlined in Section [IV Al 
we ignore the derivative operator in V2 and calculate the 
correlators of the corresponding primary operators. The 
pertinent correlation function decomposes as 

(J] Vi{zi,Zj) J] V 2 {z a ,z a )) = (50) 
jeii aei-2 

(TT e l ^k Vl(Z] ' Sj) Yl e^^-'^XTT e l ^ V2{Za ' Za) ). 

jeii aeli aeli 

The correlator involving (p x contributes the following 
holomorphic factor 

II M^\rf n M*ab\-r)* 

i<j£li a<beh 

x n M^alrfg^A^), (51) 

iS/l ,06/2 
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where Ql ts {z) = ?i^} s (z). This notation is chosen to 

be consistent with the general hierarchy, where 7i^\ s is 
used for the correct CM dependence at level n. Whenever 
convenient, we will suppress the r and t dependence and 
drop the bar on s, still keeping track of the proper range 
of this index. 

The second correlator contributes a factor, which can 
be parametrized as 



I] tfiMr)i£ s 2 (5Z 2 ), 



(52) 



a<bel 2 



with 



G 2 S (5Z 2 ) = ^(_l)*fe e ^15(fe+^) 2 e 2 7 r 4 (fe+a 2 )5Z 2 (53) 
fc 



and a 2 



15 



The basis set of wave functions is 



then obtained as a product of the conformal blocks ([5"Tj) 
and (|52p . The derivatives commute with the magnetic 
translation operators and (jTTJ) must be fulfilled for all 
divisions into sets separately. Therefore, we can consider 
the functions 



(54) 



instead of the full wave function without loss of general- 
ity. The CM part, H {2) (Z^\ is a sum of products 
gl x (ZW)g* 2 (ZW) and ip jas denotes the Jastrow factor 
expected from the result on the plane. Even though 
Z^) depends on the CM of the sets and, 
therefore, on both total CM and relative coordinates, we 
will still call it the CM function in the following sections. 

Note that we choose to express the correlators in the 
basis defined by (|42|) rather than the one obtained di- 
rectly form p6p . This will be very convenient, since by 
construction the basis we use already incorporates the 
correct boundary conditions for the first set. As we will 
see, a similar thing happens for a general hierarchy state, 
where a proper choice of conformal blocks will automati- 
cally impose the correct transformation properties for all 
but the last set of particles. 

As for the Laughlin states, we proceed by demanding 
t a \ to be diagonal (with eigenvalues +1) on the wave 
function. As the transformation properties of both the 
Jastrow factors and the gaussian are known, we focus 
only on the CM dependence and derive its properties 
under single-particle translations of %i <E I a - 



nv>{z a ) 



(55) 



We suppress the dependence on all CM coordinates but 
the translated one in this and the following section. We 
also introduce |Q(")| 2 abbreviation. Note, 

that K a is odd for all a's, since the vertex operators are 
fermionic. Thus, all electrons obey the same boundary 
conditions, independent of which set they are in. 



bmce is a linear combination of products Q\ Q , 
we first determine their properties under single-particle 
translations. As already pointed out, the former has, 
by construction, the correct properties for translations of 
particles in the first set, but not for those in the second: 



2-jriai 2 



gjl + ClT ) - p -2™(2Z 1 + fZ 2 ) - J7 rr§ /R fi x 

Q] 1+2 {ZW) ~ ■ (5bj 

Comparing these expressions with the ones of Q^ 2 (Z^): 



= 2iria 2 5 



Gl(Z^+q2) 

Gl 2+5 (Z^) ' [ ] 

we find that the coordinate dependent factors (indepen- 
dent of the free parameters s\ and s 2 ) combine to give 
2Zi + 3Z 2 = Q( 2 ) • J2 a Q (q) Z a and, thus, the correct 
factor to cancel the one from the relative part and the 
gaussian. This is not a coincidence, as it may seem here, 
but a general property due to the construction of the ver- 
tex operators. Correct transformation properties along 
L\ require special combinations of s\ and s 2 , so that 
(2si + s 2 )/3 is an integer. Hence, we find that only the 
combinations where —s± + s 2 = mod 3 are consistent 
with the eigenvalues in (|55[) . This reduces the number of 
'good' basis states from 45 to 15 (x 4 for the different 
flux sectors). By taking r = (N s — «i) mod 2 we obtain 
the correct sign in ([55| . 

We can now use this reduced set to find the linear 
combinations that also transform correctly under t 2 l . It 
is easy to see that if a pair (si, s 2 ) satisfies —si + s 2 = 3k, 
then (si + 2, s 2 + 5) does, too. Therefore, t 2 ° only maps 
functions in this reduced set into each other. Hence, we 
use the parametriztion (si,s 2 ) — (21,51 + 3s') and form 
the linear combination 



1=0 

x gI l+ 3.AZ {2) ), (58) 



which transforms correctly, if t is chosen to be: t — N s — 
«i. Note, that we use s = 3s' mod 5, s' — 0, ... ,4 to 
label the CM coordinate. That this is a natural choice is 
seen by the action of T\ on ([551) . 



2 



_ ifg^)g 2 l+3 Az^) 



1=0 



x e 27rm 1 (3+2)2$- e 2™a 2 5 2$ 7 

= H<?\zW,zW)e 2 ™i% 



(59) 



where s naturally occurs multiplied with the filling frac- 
tion v — 2/5, which amounts to K\ — (2 As) mod N 2 . 
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(2) 

In addition, T2 acting on 7i g only shifts s — > s + 1. 
Thus, H.^ is invariant under Xf as expected and .K2 = 0. 
There is no overall sign, as both Q 1 and Q 2 pick up a fac- 
tor (—1)' under s — > s + 5. 

In order to get the wave function, we need to reintro- 
duce the derivatives and antisymmetrize over all possible 
divisions into the two sets: 

*2 /5 = e (-i^n^ n m^jt) 3 

»1<»2<--»iv/2 fc ij<il 

a l <a 2<- a jv/2 

X LI MZa jai \T) 3 Y[M^ ai \T) 2 

x Hf\z^\Z^)e-^^^ (60) 

with Z« = S'EjZii/Li + I^ZajLx and = 
5j^-;2 aj /Li. There are \N s /2 — 1] possible choices for 
the derivative operator for a fixed set of quantum num- 
bers of ^2/ 5- More comments on this can be found in 
Section I VII where we test the candidate wave function 
numerically against exact diagonalization results. 



C. The general hierarchy state 

In this section we generalize the previous construction 
of the Laughlin and v — 2/5 wave functions on the torus 
to the general class of filling fractions discussed in Section 

EH 

Recall that to construct the hierarchy states at level n 
the electrons are divided into n sets I a . The electrons in 
set I a are represented by Vq(oo , giving a center of mass 
coordinate Z a — J2kei a Z k/Li- Thus, the correlator we 
want to compute is of the form: 

n 

(IlllV'fe^- ( 6l ) 

a=l kel a 

It factorizes into a product of correlators containing only 
one field, each of which can be calculated as in Section 
IIIII In the following, we use the same conventions as in 
the previous section. We denote the functions Grts as 
CM function, even though they depend on the CM of the 
sets and not the total CM. For simplicity, the fluxes are 
set to zero, 4>i — 4>2 — 0, and we suppress the indices r 
and t as they are fixed by the boundary conditions. 

The wave functions are extracted from (|6ip by ex- 
actly the same approach as used in the Laughlin case 
and v = 2/5. As in the latter, the wave function is 
an antisymmetric sum over the various ways to divide 
the electrons into the sets. Each summand must be an 
eigenfunction of the magnetic translation operators. This 
constraint is used to construct the suitable linear com- 
binations of the conformal blocks for any given division. 



The basis states are given by 

n 

*r.«.{..}=n n m-w» q(q, - q(q) n 

a—l i a <j a el a a<f3 
n 

a=l 

where we used the same conventions for Q J S . as in the 2/5 
case: 

fc 

with a j ■ — + ^ and Sj = 1, . . . AR 2 — 1. Their transfor- 
mation under translations along L\ and L 2 follow directly 
from (155)) : 

— —. — = e j e j , (64) 

and they have the periodicity property 

Gi. +R .{Z$) = {-lfe-^e-^^gi^) . (65) 

The functions Qi depend only on Z$ = Cj J2a=j+i + 
qjZj and are thus independent of the first j — 1 CM 
coordinates. In addition, they are identical for the parent 
and daughter state, except that Z$ must be replaced by 

(i) 

Zn+i- This yields a recursive construction of the full CM 

dependence Ti-i^ . To be more specific, we will show by 
recursion that at filling fraction v n = p n /q n there are 

exactly q n possible CM functions Tig given by: 

h£\zW,...,zW) (66) 

<Ji — 1 qn-l— 1 n 

= £■■• E II (-i)^ + - + ^gi(z^), 

z 2 =o ;„=o a=l 

where 

n 

Sa = qJa + E ° alk ' a < 71 
k=a+l 

s n = qJn + q n -is . (67) 

In analogy with the v = 2/5 case, and for reasons that 
will become clear below, we choose s n as the free param- 
eter, with s n = q n -is mod q n , s = 0, . . . , q n — 1. 

We have already shown that the CM part of the Laugh- 
lin wave functions is of the form (|66p . It remains to 
show that, assuming (|66[) for a parent state at level n, 
it follows that also the CM part of the daughter state 
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at level n + 1 obeys this relation. To go from level n 
to level n + 1 in the hierarchy, we consider a quasiparti- 
cle density of l/(2fcj + i + 1) atop the parent state. The 
daughter state is then at filling fraction z/ n+ i determined 
by (flC)]) . To construct the wave function, there is then 
an additional vertex operator, V n +i, with charge vector 
Q(n+i) = ... 5 jn.^ g_L). Note, that all charge vec- 
tors are now n+1 dimensional objects. The wave function 
will be of the form 

where Ti^t^ is some suitable linear combination of prod- 
ucts of Qg . 

Invariance under one-particle magnetic translations 
puts constraints on the wave function. As the Jas- 
trow type factor is known, we can compute which re- 
lations H.g™^ must satisfy under one-particle (normal) 
translations. Under the translations Z a — > Z a + 1 and 
Z a — > Z a + t, Wj™^ transforms as 



x e 



(_l\(N s -K a ) 

(69) 



where we use the same simplified notation as in (|55|) : 

+ C ) = ^(*. Z n + W but with Z a 
replaced by Z a + c in the Z^'s. In deriving the second of 
these relations, we used (fT4|) . Equation ([69]) is valid for 
the CM dependent functions at all levels, in particular 
also at level n. 

The function H^t^ is a linear combination of prod- 
ucts: n"=i ^f(^n+i)- Thc task of finding the correct 
linear combination is simplified by using the result for 

the parent state. By construction, (Z^+j, — , ^n+i) 
transforms correctly under translations of all CM coordi- 
nates Z a but Z n+ i . For this last coordinate, note that 
for a < n, G^ depends in the same way on Z n+ \ as on 
Z n , as = q a Z a + c a (Z a+ i + ... + Z n + Z n+ i). This is 

not true for Gsn , but we can separate the difference by 
writing: 

SjZHzQj = g£>(q n (Z n + Z n+1 ))f(Z n ,Z n+1 )(70) 
with 

tfr? 7 \ _ Ss n (QnZ n + C n Z n+1 ) 
f(Zi n ,Z, n+ i — — — ■ -. (II 

y% n (q n Z n + q n Z n+ i) 

This implies that for a given set (l 2 , ■ ■ ■ , l n ) in (|6"6"|) . the 
product Q 1 ■ ... ■ Q" transforms in the same way under 
translations of Z n and Z n+ \, except for the ratio (|7ip . 



The transformations of this ratio is, however, readily ob- 
tained from (|6"5|) , and since it is independent of l n , we can 
infer the following relation for translations along L%: 



n { ^(z n+1 + i) 



^nt\z n+1 ), (72) 



and a slightly more complicated relation for translations 
along Li: 

H { ::\z n+1 +r) 

-2™E? =1 (Q ( " + 1) -Q tfJ % + J-^ + i) -i Trr ( K „ + %2i) 

= e 3 e B " 

qi — l 971-1 — 1 n 
l 2 =Q ;„=0 a=l 



.too 



= ±e i 3 e R " 



xH { Z } +c (Z n+1 ). 



(73) 



The sign is calculable, but is not needed in the following 
discussion. As q n and q n +i are relatively prime, jc n mod 
q n runs through all possible values of q n . Therefore, in- 
variance under translations of the set I n +i will require 
linear combinations of all q n CM functions of level n. 

To obtain an eigenfunction of , Ti.^ must be com- 
bined with C^i+i ) w fii°fi transforms according to 
([64]) with c = q n+ i : 



y S „ + l\ Z n+l +Qn+1) 

pn+l / 7 (n+l)\ 
y 8n+1 {6 n+ i ) 



2iri- 



(74) 



— — ~, ' rr — = e »+i e 

Comparing the first lines of (|69|) . ([72]) and ([74]) we see 
that the condition 



- s n + s n+1 = mq n 



m integer 



(75) 



must be fulfilled. This reduces the number of allowed 
products of conformal blocks from q n Rj l+1 to -R^ +1 . A 
further reduction is obtained from requiring invariance 
under translations along the r-direction. The only com- 
binations that transform correctly are given by 

nttl\zi%,...,Z^) (76) 
g„=i 

= (-l) tj n^(Z n+1 + 3 r)G n q +s{Z n+1 + jr), 
j=o 

where, again, s n +i = q-nS and s = 0, . . . , q n +i — 1- Insert- 
ing ([75]) into the above equation, it is straightforward to 
verify that it is of the form (f6"6")> . We chose s n = in (jT6"|) 
in accordance with the sign convention in (|66|) . Taking 
another value amounts only to an overall sign change. 

We can also compute the quantum numbers of the wave 
functions (|68p under magnetic translations recursively. 
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Under translations with T\ , the conformal blocks at level 
n + 1 pick up a phase given by: (_i) r ("+ 1 ) e 27rw ™+i ;s >«+i . 
Thus, K\ is given by Ki = (r(n + l)JV a /2 +p n+ iNs n+1 ) 
mod AT S . We can also show that K 2 = t(n+l)N s /2 mod 
N s . As in the previous cases, the different fluxes cj)^ in 
(|17|) can be incorporated by a proper choice of r and i. 

Hence, we found exactly q n +i conformal blocks with 
the correct quantum numbers. The reader may also note 
here that the total number of conformal blocks (before 
imposing the boundary conditions) is not given by the 
filling factor. Other choices for the charge vectors may 
decrease this number considerably and simplify compu- 
tations. This will be discussed in more detail in the next 
section. 

In conclusion, we have constructed the hierarchy wave 
functions in torus geometry for all filling fractions v n 
that are obtained by successive condensation of quasi- 
electrons. At level n we find q n wave functions that are 
eigenf unctions of the magnetic translation operators Xi 
and T 2 n , confirming the expected ground state degener- 
acy on the torus. Reinstalling the derivatives and an- 
tisymmetrizing, we obtain an explicit expression for the 
wave functions at filling factor v n : 



= e (-D-n n^n n 

h,...,I n k=2i k el k a<0 

n 



1 ,, 2 



(77) 



The sum runs over all possible ways to divide N parti- 
cles into n sets of size Mi, . . . , M n and (— l) a is the sign 
picked up by rearranging the radially ordered particles 
into the sets. 



V. ALTERNATIVE CHARGE VECTORS 

We have already mentioned that the charge vectors 
Q( Q ) are not uniquely determined by the wave function. 
To spell out this ambiguity we first notice from (fr?2")) that 
the charge vectors enter in the relative part only through 
the scalar products Q'"- 1 ■ If we introduce the vec- 

tor c a = 1/R a , the holomorphic current operator takes 
the form: J{z) = i Y) a c a d<p a {z). The constraints on the 
vertex operators that are imposed by the short-distance 
behaviour of the electrons, unit U(l) electric charge and 
the filling fraction, can be expressed as the following con- 
ditions on the n-dimensional vectors Q^ Q ' and c, 

q(«) . q(/3) = K tP 

Q( Q ) • c = 1 

c-c = v. (78) 

These scalar products are invariant under simultaneous 
0(n) transformations on and c. (It should also 



be clear how to include quasiholes by introducing a new 
set of vectors, l^ a \) We note, however, that an 0(n) 
transformation will in general give irrational radii imply- 
ing that the CFT is not rational. This means that the 
charge sums cannot in any obvious way be reorganized 
as a finite sum over conformal blocks as required by the 
general analysis in Section IIIII Naively it would seem 
that no similar restrictions to rational radii would apply 
when working in the plane, but it is not unlikely that 
they would emerge from a study of the edge theory for a 
finite droplet. 

We have not tried to find the most general transforma- 
tion that leaves (|78[) invariant and still maintains rational 
radii of the fields, but only studied a few examples. Note 
that it is not obvious that preserving (|78[) automatically 
implies that the CM part of the wave function remains 
the same. On the contrary, this amounts to non-trivial 
mathematical relations between conformal blocks on ra- 
tional tori. We have numerically verified such a relation 
explicitly in the case of v = 2/5 and v = 3/7, and suggest 
that similar results hold in general. If, on the contrary, 
these relations turned out not to hold, it would indicate 
the presence of inequivalent hierarchy states at the same 
level, and with the same filling fraction. 

We now present alternative charge vectors for the 
fractions 2/5, 4/11 and 3/7, in a basis where c = 
(v 7 ^, 0, . . . , 0). In this basis a background charge is 
needed only for the field tpi, since the correlators for the 
remaining fields are neutral by construction. As we shall 
see, the number of conformal blocks is much smaller in 
this basis, which simplifies calculations. On the other 
hand, the hierarchy structure is not manifest, e.g. the 
computation of v = 2/5 blocks is in no simple way re- 
lated to the computation at v = 1/3. 

For v — 2/5, the charge vectors in this basis are given 
by the following symmetric expressions, 



n(i) = (A L) 



Q (2) = ( 



10 



(79) 



where the sets I\ and I 2 each contain A^/2 electrons. 
Since the radii squared are even integers, the double 
charge-lattice sums reduce to single sums (p- = 1), and 
the conformal correlator gives a holomorphic basis of the 
form 

i<jell a<beh 

x I] ^i(z w |r) 2 F s 10 (5Z)F r 2 (Z 12 ), (80) 
where Z = Z\ + Z2, Z\i = Z\ — Z2, and 



F™(5Z)=J2< 



,i7i-rl0(fc+s/10) 2 2iri(k+s/10)5Z 



e v !»/■"•/ g 

Fl(Z 12 ) = ^ e i ' rT2(fc+s ' /2)2 e 27ri( ' £+s ' /2)Zl2 . (81) 
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The first conformal block depends on the total CM, while 
the second contains only relative coordinates. However, 
it is still not possible to separate the CM dependence, 
as explained below. A linearly independent basis of 20 
states is obtained by taking s — 0, ... ,9, s' = 0, 1. This 
should be compared with the calculation in Section 1VB, 
where the original 720 states had to be reduced to 20 by 
imposing proper boundary conditions. 

To compare with our previous calculation we change 
the basis for the <pi field to 

F W t (5Z) e inTlO(k+s/5+t/2) 2 e 2TTi(k+s/5+t/2)5Z ^) 

k 

where the parameters have the values s = 0, . . . , 4 and 
t — 0, 1. In this parametrization, s labels the CM coor- 
dinate (in fact, this is precisely s' in section HV B|) . while 
the boundary conditions are coded in r and the differ- 
ent combinations of t and s'. An analysis along the lines 
of that given in the previous examples, shows that the 
following linear combinations give eigenstates of t t | and 



Hi 2 \Z 1 ,Z 2 ) = Fl%{5Z)Fl{Z 12 ) 

+ (-1)PF™(5Z)F* +1 (Z 12 ), (83) 

where we choose a — (N s — ki + 4>i/tt) mod 2 and 
(3 = (N s — ki + 4*2/ k) mod 2. Note, that the boundary 
conditions require a non-trivial combination of different 
CM functions and, thus, there is no simple way to fac- 
tor out the total CM. As already mentioned, we have 
checked that the wave functions thus obtained are nu- 
merically equal to those given in (|55|) . This amounts to 
a rather complicated identity between sums of products 
of generalized theta functions. 

Finally, we also give explicit expressions for the charge 
vectors for v = 4/11 and v = 3/7. In the former, we have 
different number of particles in the two sets and, thus, 
the charge vectors look less symmetric: 

Q (1) = (/f,~), Q( 2 > = (/f,-^). (84) 

For 3/7 we find the three charge vectors: 




Using these, the number of conformal blocks is greatly 
reduced compared with the earlier representation (|12|L 
but unlike the v = 2/5 case, there remain states that 
must be excluded by applying the boundary conditions. 



VI. NUMERICAL TESTS 

We have carried out numerical tests of some of the sim- 
plest hierarchy wave functions, namely the v — 2/5 Jain 
state, and the recently discovered v = 4/11 state, both at 
level two, and the level three Jain state at v = 3/7. Here 
we present results for the 2/5 state and just comment 
briefly on the others at the end of the section. We com- 
pare our wave functions to the ground states obtained 
by exact diagonalization of eight respectively ten parti- 
cles using an unscreened Coulomb interaction. 

As discussed in Section IIII B[ the wave functions we 
have constructed are not unique because of the freedom 
associated with the derivative operators T)( a K In par- 
ticular, for the 2/5 state we have to define the operator 
X>W. A simple set of choices is, 

x£> = q (86) 

and we shall denote the corresponding wave functions by 

(n) 

ip2/5- For the simplest choice of n = 1, we find that 

the overlap between the exact groundstate and ^2/5 i s 
well above 0.9 for all values of L\ and N < 10. For some 
values of L\ this overlap even exceeds 0.99, see the lowest 
lying curves in Figs [T] and [2j 

The results are improved further if we take a linear 
combination of states ^2/5 with different values of n. We 
find that the overlap with the exact states and the space 
spanned by {"02/5 }n=i f° r different values of k quickly 
converges to a number very close to one — in fact, above 
0.99 already when two or three states are taken into 
account — irrespective of the values of L\ and TV < 10. 
These results for eight and ten electrons are summarized 
in Figs [1] and [H We note that the small variation of the 
overlap as a function of L\ very much resembles results 
obtained earlier for the Laughlin state on the torus (26| . 

We find the fact that we only need to take two or three 
states into account to get a wave function as good as the 
Laughlin state as strong evidence for the correctness of 
our approach, especially in the light of that the Hilbert 
space in the K = sector for ten particles at v = 2/5 
has 26152 dimensions. 

Preliminary results for the state at v = 4/11 — one 
of the newly observed non-Jain states — also seems very 
promising. However, both the numerical studies and the 
interpretation thereof are more involved and will be pub- 
lished elsewhere. Moreover, we have compared the trial 
wave function at filling fraction v = 3/7 with exact diag- 
onalization results and found that already the simplest 
choice of the derivative operators (|24[) gives a reasonably 
good description. 

VII. SUMMARY AND OUTLOOK 

To summarize, we have constructed torus versions of 
newly proposed wave functions describing the Haldane- 
Halperin hierarchy of incompressible quantum Hall states 
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Figure 1: Projection (amplitude) of the eight-particle exact 
solution (obtained by diagonalization of unscreened Coulomb 
interaction) to the subspace spanned by {ipys}n=i , with k = 
1 •■>■ 



example, the Halperin (m,ri2, m) states. It should also 
be possible to construct quasiparticle states using simi- 
lar techniques. For the quasiholes in the hierarchy this 
should be straightforward, since they are described by lo- 
cal vortex operators that are primary fields in the rational 
CFT's used in our construction. The quasielectrons are 
more difficult, but a suitable adaptation of the methods 
developed in 0] is likely to work. We also believe that 
the quasielectron excitations of the Moore-Read pfaffian 
state can be obtained. 

Finally we note that our construction provides a way 
of investigating the adiabatic continuity from the solv- 
able thin torus case to the bulk for the hierarchy states 
considered here. 
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Figure 2: Projection of the ten-particle exact solution to the 
subspace spanned by {^2/s}n=i> with k = 1, . . . , 12. 



in the lowest Landau level. In particular, we managed 
to incorporate the homogeneous background charge, and 
derive the non-holomorphic gaussian factor, in a mathe- 
matically sound manner. In the previous calculations in 
the disk geometry boundary terms had to be ignored @ , 
and in a spherical geometry one typically put a compen- 
sating charge at infinity at the price of not obtaining the 
gaussian factor. 

At a technical level we note that all the boson radii 
R a are square roots of integers. One can construct wave 
functions based on CFT where (some of) the i? a 's are 
rational, but these states do not seem to be part of the 
usual hierarchy, and might even be non-abelian. 

We believe that the techniques developed in this paper 
should be useful also to construct other QH states, as, for 



Appendix A: THE BACKGROUND CHARGE 
CONTRIBUTION 



In this Appendix we evaluate the integrals in (127|) and 
that originate from the homogeneous background 
charge. We consider a general rectangular integration 
domain defined by the complex number £ = a + ib 



d 2 j 



L+Li 



dx 



b+L-2 



dy. 



(Al) 



The integral in (J27J) is now easily evaluated, and, when 
combined with the second exponential yield factors, 



exp 




q a Z a - N s 



Li 



1 



(A2) 



First notice that for the choice £ = {L\ + iLi2)/2, corre- 
sponding to the symmetric integration region used in the 
main text, the integral in (|27[) vanishes. For an arbitrary 
£ the effect of the integral can be incorporated by shifting 
the coordinates Zi —* Z\ + {-f ^p~) provided 



Va : 



}^M a = C a N s . 



(A3) 



These relations, which determine the M Q 's, were derived 
previously in the plane as a consequence of the assump- 
tion of homogeneity. It is to be noted that exactly the 
same relations are needed on the torus to guarantee that 
the CM part of the wave function factorizes. Without 
starting from the plane we could thus have obtained (|A3|) 
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as a consistency condition directly on the torus. Also 
note that since (|A3|) determines the M a 's, it also im- 
plies the relation (|14ll that was crucial for obtaining ex- 
plicit expressions for the CM dependence for the hierar- 
chy states. We now turn to the integral in (|28|) 



J C CM) = / dVln 



■\r) 



(A4) 



Up to a factor independent of z, this integral can be 
calculated by using the quasiperiodicity ([57)) of the theta 
function under lattice translations: 



a+Lx 

dx' 

a-\-L\ — x 



b-\-L 2 I _ 

dy' Inltf^-— |r)| s 



z 



rb+L 2 -y 

dx' / dy' ln\M^\r)\ 2 

a+Li pb-\-L 2 

dx' 



z 



dy'ln|^(^|r)| 2 



(A5) 



f dx' [ dy' ( 4ny- + 27rlmr 

Ja Jb-v \ Lx 



k(0) + 2n 



n 2 



[y-(b 



L 2 f 



The third identity follows from that the integrand is in- 
variant under lattice translations along the x-axis, while 
under z — > z + r it picks up an additional term, 



In \ e -^r e -2^z/ Ll |2 = 47ry / il + 27r l mr . 



(A6) 



Again taking £ = — (£i + iL 2 )/2 we recover (|30| in the 
text, and taking an arbitrary integration domain, i.e. an 
arbitrary £ just amounts to shifting the coordinate sys- 
tem. 



Appendix B: DERIVATION OF EQUATION (j33[) 

We assume a compactification radius i? 2 = 2p/p', 
p,p' € Z — only in this case do the vertex operators 
(with integer charges) define a rational CFT — and con- 
sider the following sum: 



i7rra _ — i7rra m 2-Ki \a e 



E 

e,m 



Thus, we can rewrite the original sum (|B1[) as 



^ ' ^ ' e iTTT(2pp' kx+sp' +rp) 2 /2pp' g ^(2pp' fci+sp' +rp)Z (a) 



r.s fei 



E 



e -i7rr(2pp'fe 2 +s!3'-i , p) 2 /2pp' e - i ^(2j9p'fc2+sp'-' , p)Z (a) 



or, introducing additional notation, as 

p'-l 2p-l 
r=0 s=0 

where 

^ S (^ (Q) ) 

_ 'y ' e i7rr(2pp'fc+sp'+rp) 2 /2pp' e ^(2pp'fc+sp'+rp)Z (<l) 



(B3) 



(B4) 



(B5) 



This is (|33|) . note that Eq. IB4l is a finite sum of factorizcd 
terms. 



Appendix C: THE v = 3/7 AND 4/11 STATES 

As a service to the reader that does not want to pen- 
etrate the general formalism of Section IIV1 we here pro- 
vide explicit expressions for the v = 3/7 and 4/11 wave 
functions. 

The latter is, as v = 2/5, a level two filling fraction, 
but with a quasielectron density of 1/3. The computation 
of the torus wave function is in complete analogy to 2/5 
when is replaced by: Q^ 2 ' = (-^=, -±L). The fact 
that now the sets M a have different size is of no conse- 
quence for the calculation. We find Mx = 3M 2 = 3M for 
N = AM electrons. Using ([53]), but with radius R\ = 33 
instead, (|58|) takes now the form: 



n,(2) (z (l) 7 {2) s _ V^/ ,stZ/jl ( 7 a) z (2) \ 
rtg l Z 4/ll' Z 4/lli - 2^~ 1 '' 1^4/11 ' X 4/l J 



!=0 



x Sxxi+3s'(^4/\x) 



(CI) 



The integers n, r 2 and ti, £2 are fixed by ([TT)) . For 

4(e+mJ? 2 /2)z<' l »-(e-mi? 2 /2)z( £ '»] 0i = 02 = 0, we find n = r 2 = N s - q and also tt=t 2 = 

N s — q. The total wave function can then be written as 



(Bl) 



M 



We write 

mi? 2 /2 = mp/p' = (fhp' + r)p/p' = rhp + rpjp' , (B2) 

where r = 0, . . . ,p' — 1. Additionally, we write e + mp = 
2pfci + s and e — mp = 2pk 2 + s, with s = 0, . . . , 2p — 1. 



4/11 



i l< i 2<--- i 3Af fe=l ij<h 

a 1 <a 2 <...a M 

x (^ n ,^ 2 / n )e-^^^, (C2) 
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where 



with the complete wave function given by: 



3M 



M 



^ 2 /ii = nEV^ 



(C3) 



As was the case for v = 2/5, the derivatives in (|C2[) are 
not unique and each choice yields a trial wave function. 

The filling fraction v = 3/7 is constructed by 3 vertex 
operators with charge vectors Q^ 1 ) = (-^,0,0), Q^ 2 ^ = 

^,0) and QW = ^). The sets are of 

the same size M, for N = 3M electrons. Again, we set 
both solenoid fluxes to zero, thus boundary conditions 
require t — r — N s — 3. Then, (|6"6")) takes the form: 



M M 



3/7 



e (-irn^m^? n^i- 



i l< i 2<-'M 
i 1 <a. 2 <...a M 



fe=l Z = l 



ij <ii 



a,j<ai a.j<on ij ,&i 

JlM^Mf n MZa jai \T) 2 



nj(2) (7 {l) 7 (2) 7 (3)n 
x «s ^3/7' Z 3/7' Z 3/7^ 



(C 



oy(3)/ z (l) 7 (2) 7 (3)x _ V^/ 1 Nt(i 2 +i 3 )r;i r7 (1)N i 



i 2 =0i 3 =0 



2l 2 +2l 3 \ ZJ 3/7J 
(3)- 



As explained earlier, (— l) a is the sign picked up by rear- 
ranging the ordered fermions into sets. Antisymmetriza- 
tion requires the derivatives T>^ and to be different. 
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